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ABSTRACT 

Aims. We probe the physical conditions in high-redshift damped Lyman-fy systems (DLAs) using the observed molecular fraction and 
the rotational excitation of molecular hydrogen. 

Methods. We search for Lyman- and Werner-band absorption lines of molecular hydrogen in the VLT7UVES spectra of background 
QSOs at the redshift of known DLAs. 

Results. We report two new detections of molecular hydrogen in the systems at z a bs = 2.402 and 1.989 toward, respectively, 
HE 0027- 1836 and HE 23 18- 1107, discovered in the course of the Hamburg-ESO DLA survey. We also present a detailed anal- 
ysis of our recent H 2 detection toward Q 2343+ 125. All three systems have low molecular fractions, log/ < -4, with / = 
2N(H 2 )/(2JV(H 2 ) + N(Hi)). Only one such H, system was known previously. Two of them (toward Q 2343+ 125 and HE 2318-1 107) 
have high-metallicities, [X/H] > -1, whereas the DLA toward HE 0027-1836 is the system with the lowest metallicity ([Zn/H] = 
-1.63) among known H 2 -bearing DLAs. The depletion patterns for Si, S, Ti, Cr, Mn, Fe and Ni in the three systems are found to be 
very similar to what is observed in diffuse gas of the Galactic halo. Molecular hydrogen absorption from rotational levels up to J = 5 
is observed in a single well-defined component toward HE 0027-1836. We show that the width (Doppler parameter) of the H 2 lines 
increases with increasing J and that the kinetic energy derived from the Doppler parameter is linearly dependent on the relative energy 
of the rotational levels. There is however no velocity shift between lines from different rotational levels. The excitation temperature 
is found to be 90 K for J = to J = 2 and ~500 K for higher J levels. Single isothermal PDR models fail to reproduce the observed 
rotational excitations. A two-component model is needed: one component of low density (~50 cm 4 ) with weak illumination (x = 1) 
to explain the J < 2 rotational levels and another of high density (~500 cm 4 ) with strong illumination (x = 30) for J > 3 levels. 
However, the juxtaposition of these two PDR components may be ad-hoc and the multicomponent structure could result either from 
turbulent dissipation or C-shocks. 

Keywords, galaxies: ISM - quasars: absorption lines - quasars: individuals: HE 0027- 1836, HE 23 18- 1107, Q 2343+125 

1. Introduction high-redshift (z a bs > 1.8) H2-bearing DLAs are rare with only 

ten detections report e d up to now, wi t h reds hifts up to z a bs - 4.2 
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H , High-redshift damped Lyman-a systems (DLAs) are the ab- dLedoux et al.ll2006h . iPetitiean et alJ (120061) noted* a correlation 

sorbers of highest H t!| column density seen in QSO spectra, with between me presence of molecular hydrogen and the metallicity 

N(Hi) > 2 x 10 2U cm" 2 . These absorbers are a major reser- ([X/H] = i og JV(X)/JV(H) - log X/H , with X = Zn, S or Si) in 

voh_of neutral hydrogen at 2 < z < 3 (e.g. | Prochaska et al. | the gas that may help select the best candidates for future sur- 

|2005l) and have lo ng been identified as the precursors of present veys Indee d, 40% of the [X/H] > -1 .3 DLAs have log / > -4, 

day galaxies (see |Wolfe et al.||2005[ for a recent review). Our with f _ 2jV(H 2 )/(2AT(H 2 ) + N(Ui)). On the other hand, only 

understanding of DLAs is primarily based on absorption-line 15% of the overall sample (with _ 2 . 5 < [X /H] < -0.3) have 

studies, involving the detection of low-ionization metal transi- log y > _ 4 Furthermore, up to now, there is no H 2 detection 

tions and, in a few cases, molecular hydrogen (e.g.|Ledoux et al.| in [X /H] < -1.5 DLAs, down to a detection limit of typically 

l2003tlSrianand et al.ll2005l) . Molecular hydrogen is conspicuous jV(H 2 ) = 2 x 10 14 cirT 2 
in our Galaxy: l ines of sight with log/Y(Hi) > 21 usuall y have 
log/V(H 2 ) > 19 dSavage et al.|[~1 977; Rac hford et al.ll2002l) ' whilst 

The distribution of molecular fractions is observed to be bi- 

— 7 modal: all the upper limits are lower than log / = -4.5, whereas 

Based on observations carried out at the European Southern j . .• u • ,. •. , , /> i i i i 

Observatory (ESO) under progs. ID 67.A-0022, 69.A 0204, 072.A- f detections above this limit have log/ > -3. Low molecular 

0346, 072. A-0442, 073. A-007 1 and 074. A-0201 with the UVES spec- fractions (with log/ <: -4) are basically unexplored (see how- 

trograph installed at the Very Large Telescope (VLT) Kueyen UT2 on ever |Levshakov et al.||2002|). Up to now, only two low-/ sys- 

Cerro Paranal, Chile tems have been reported, i n the DLA at Zabs = 3.025 toward 

1 To avoid confusion, we use the convenient spectroscopic notation Q 034 7-383 (log / = -5.9;|Levshakov et al.|200l|Ledoux etjp 



(e.g. Sin for Si + ) for atoms and ions following the convention widely 120031) and in the DLA at Zabs = 2.431 toward Q2343+125 
used in this field. (log f = -6.4; IPetitiean et al.ll2006l) . 
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Molecular hydrogen is primarily formed on the sur- 
face of dust grains in the interstellar medium (e.g. 
iHollenbach & Salpeterl 1 197 lb . but can also result from the 
formation of negative hydrogen (H ) if the gas is warm and 
dust-free (e.g. Bla ck et al]|l987l) . H2 is also formed in the gas 
phase by radiative association of H + with H to form H^, fol- 
lowed by charge transfer with H. However, t he reaction is slow 
in the conditions prevailing in the DLA gas. ISpitzer & Jenkins! 
( 1 1975b first showed that the observed relative populations of the 
H2 rotational levels cannot be described by a single excitation 
temperature. Whilst the low-J levels can be thermalised, the 
high-J levels are populat ed by direct formation on these levels 
and UV pumping (e.g. lJural 1 1975b or by collisions if part 
of the gas is heated to temperatures above a few hundred 
K by, for example, turbulent dissipation and/o r shocks (e.g. 
Uoulain et alj[l998l: ICecchi-Pestellini et al.ll2005h . In DLAs, H 2 
excitation is generally ex plained by predominant UV pumping 
dHirashita & Ferrarall2005l) . 

Here we present the analysis of H2 absorptions in three 
DLAs, at z abs = 2.431, 1.989, and 2.402 toward, respectively, 
Q2343+125, HE2318-1107, and HE0027-1836. The three 
systems have low molecular fractions according to the above cri- 
terion, log / < -4. The z abs = 2.43 1 system toward Q 2343+125 
has the lowest molecular fraction measured till now in a DLA. 
We recently reported this detection as part of a high-met allicity 
([X/H] > - 1 .3) UVES DLA sample dPetitiean et al.l2006l) . A de- 
tailed analysis of the system is given here. The latter two systems 
are ne w detections discov ered during the Hamburg-ESO DLA 
survey (Smette et al. 2005, Smette et al. 2007, in prep.) and are 
reported for the first time. H2 is detected in rotational levels up 
to J = 5 in the DLA toward HE 0027-1836. This is the system 
with the lowest metallicity amongst known H2-bearing DLAs. 
We present the observations in Sect. 2 and properties of the three 
systems toward respectively Q 2343+125, HE 2318-1 107 and 
HE0027-1836 in Sects. 3 to 5. We discuss the H2 excitation 
toward HE 0027- 1836 in Sect. 6 and comment on metallicity 
and depletion in Sect. 7. We conclude in Sect. 8. 



2. Observations 

The quasars Q 2343+125, HE 2318-1 107, and HE 0027-1836 
were observed in visitor mode wi th the Ultraviolet a nd Visible 
Echelle Spectrograph (UVES, iDekker et all |2000) mounted 
on the ESO Kueyen VLT-UT2 8.2 m telescope on Cerro 
Paranal in Chile. Q 2343+ 125 was observed on October 29, 
2003 under program ID. 072.A-0346. HE 2318-1 107 and 
HE0027-1836 were observed on September 16, 17, and 
18, 2004, under program ID. 073.A-0071. HE0027-1836 
was re-observed on October 8 and 9, 2004, under program 
ID. 074.A-0201. Observations were supplemented with UVES 
data from the ESO archive, progs. 67.A-0022 and 69.A-0204 
(P.I: D'Odorico) for Q 2343+125 and prog. 072.A-0442 (P.I: 
Lopez) for HE 2318-1 107 and HE 0027- 1836. The d ata were 
reduced using the UVES pipeline (Ball ester et alj2000l) which is 
available as a context of the ESO MIDAS data reduction system. 
The main characteristics of the pipeline are to perform a precise 
inter-order background subtraction, especially for master flat- 
fields, and to allow for an optimal extraction of the object sig- 
nal rejecting cosmic rays and performing sky subtraction at the 
same time. The pipeline products were checked step by step. The 
wavelength scale of the reduced spectrum was then converted to 
vacuum-heliocentric values and the portions of the spectrum cor- 
responding to different settings were each rebinned to a constant 
wavelength step. No further rebinning was performed during the 



analysis of the whole spectrum. Individual 1-D exposures were 
scaled, weighted and combined together. Standard Voigt-profile 
fitting methods were used for the analysis to determine column 
densi ties using the va cuum wavelengths and oscillator strengths 
from Morton! d2003b for metal ions, except for the oscillator 
strengths of Fe O/U/U 1 063, 1064,1 1 12,1 121 dHowk et alJ feOOO) 
and Ni 0^1317,1370 dJenkins & Tripp! 120 06), and the wave- 
lengths of So ^1250,1253 (|Morton|[l99!|). We used the oscil- 
lator strengt hs from the Meudon groupo based on calculations 
described in lAbgrall et alj d 1994b for H9. We adopted the Solar 
abundanc es given b y Morton! (|2003b based on the meteoritic 
data from iGrevesse & Sauvall d2002l) . For each system, the ori- 
gin (v = kms~') of the velocity scale is set to the redshift of 
the H2 component. 

3. The system at z abs = 2.431 toward Q 2343+1 25 



This 
then 



DLA has firs t been st udied by ISargent et al] dl988b 



bv lLuetalJ _ J1996h^Jp;Odoricoetal.l d2002i) and 



iDessauges-Zavadskv et al.l d2004b with high spectral reso- 
lution (~ 6.5 kms 1 ). We recentl y reported the detec tion of 
molecular hydrogen in this system (Petitiean et al. 2006). From 
Voigt-profile fitting to the H 1 Ly-a, [5 and y lines, we find that 
the damped Lyman-a component is centered at z a bs = 2.431 
and the column density is log N(H 1) = 20.40 ± 0.07, consis- 
tent with previous measurement by iD'Odorico et al.l {2002): 
logiVTHi) = 20.35 ± 0.05. The fit to the damped Ly-a line is 
shown in Fig.Q] 
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Fig. 1. Fit to the damped Lyman-o' line at z a b s = 2.431 toward 
Q 2343+125. The derived neutral hydrogen column density is 
logAf(Hi) = 20.40 + 0.07. The dashed curves correspond to the 
best-fitted Voigt profile (black) and the associated errors (red). 



3. 1. Molecular hydrogen and neutral carbon content 

The optically thin H2 absorption lines at z a bs = 2.43128 are 
very weak, i.e. close to but above the 3 <x detection limit (see 
Fig. |2). A very careful normalization of the spectrum has been 
performed, adjusting the continuum locally while fitting the line 
profiles. We detect the first two rotational levels in this system 
(see Fig. 0. The total H2 column density, integrated over the 
J = and J = 1 rotational levels, is measured to be log N(H2) = 
13.69 + 0.09 (12.97 + 0.04 and 13.60 + 0.10 for J = and J = 1 
respectively). This leads to an extremely small molecular frac- 
tion log/ = log2jV(H 2 )/(2jV(H 2 )+jV(Hi)) = -6.41 + 0.16. We 
also derive a 3 <x upper limit on the detection of the J = 2 level, 
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log/V(H 2 , J=2) < 13.1, see Table Q] Neutral carbon is not de- 
tected resulting in an upper limit at 3 o~. log 7V(C i) < 12.1. This 
leads to a ratio logAf(Ci)/7V(H 2 ) < -1.6. 

Table 1. H 2 and Ci column densities in the z a b s = 2.431 DLA 
toward Q 2343+ 125. 
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Zabs 


level 


logjV(H 2 , J)* 


V 

[kms" 1 ] 


T d 
1 0-J 

[K] 


3 


2.43128(3) 


J = 
J = 1 
J = 2 


12.97 ± 0.04 
13.60 ±0.10 
< 13.10" 


1.3-2.6 


99Q+173 

< 420 




# 


Zabs 


logA'(Ci) 


logAT(Cf) 


b 

[kms -1 ] 


Avci/H, 

[kms"'] 




2.431 


< 12.10* 


< 12.40 8 
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Fig. 4. Overall Si iv and C iv absorption profiles in the z a b s = 
2.431 DLA toward Q 2343 +125. The vertical lines mark the red- 
shift of the main low-ionization metal components (see Table|2]i. 
Note that the Civ and Siiv absorptions from the H2-bearing 
component (at v = km s~') are absent (or weak). 



8 The numbering refers to that of the metal components (see Table|2]l. 

Errors represent the allowed range for logA^ITJ), not the rms error 
from the fit. 

c Depending on the inital guess, the fit converges to two values of b, we 
therefore give the corresponding allowed range. 

d The excitation temperatures between rotational levels and J are cal- 
culated using the allowed range for _/V(H 2 ,J). 
e 3 a upper limit. 
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Fig. 3. Velocity plots of observed H 2 absorption lines from the 
J = and J = 1 rotational levels at z a b s = 2.43128 toward 
Q 2343+125. The model fit is over-plotted. 



its red side. We assume that most of the metals are in this clump 
as it is associated with by far the strongest Oi absorption fea- 
ture. This main clump is itself decomposed into five components, 
among which the third one from the blue, at z a b s = 2.43129, 
shows associated H 2 absorption. About 35% of Si n and S n re- 
side in this component. The Voigt-profile fitting is shown on 
Fig. [5] and the corresponding parameters are given in Table [2] 
Note that there is very little C iv and Si iv absorption at the veloc- 
ity where low-ionization metal absorptions are the strongest (see 
Fig. |4j. High-ionization species are seen mostly blueshifted up 
to 300 km s _1 from the main low-ionization clump. It is apparent 
that the profile is much smoother far away from this component, 
suggesting a hot gas outflow. 

4. The system at z abs - 1.989 toward HE 231 8-1 107 

From Voigt-profile fitting to the Hi Ly-a and Ly-/3 lines, 
we derive a total neutral hydrogen column density of 
log/Y(Hi) = 20.68 ± 0.05. The fit is shown on Fig.0 

4.1. Molecular content 

Our spectrum of HE 23 18- 1107 covers only a few H 2 lines, 
but these are unambiguously detected in a single component at 
Zabs = 1-98888 with absorption lines from the rotational levels 
J = to J = 2 (see Fig. |7}. Some lines from the J = 3 level 
may be present but they are always blended or in regions of low 
SNR. We therefore derived only an upper limit at 3 <x: log A^(H 2 , 
J=3) < 14.7. The results of the fit are given in Table [3] The 
lines are optically thin and the quality of the spectrum is high 
enough to measure accurate column densities in each of the de- 
tected levels. The normalization of the R ~ 54000 spectrum has 
been done locally around each line. The total H 2 column density 
is log 7V(H 2 ) = 15.49 ± 0.03, leading to a low molecular fraction 
log/ = -4.89 ± 0.08. Because J = 1 is mainly populated by re- 
active collisions, the kinetic temper ature should be c lose to the 
excitation temperature 7oi - 188 K dRov et al.l l2006). 



3.2. Metal content 

The low-ionization metal profile is spread over more than 
250 kms -1 (Fig.HJi but shows a strong low-ionization clump on 



4.2. Neutral carbon 

Neutral carbon is detected in two fine-structure levels of the 
ground state, 2s 2 2p 2 3 Po (C i) and 2s 2 2p 2 3 Pi (C i*), at the same 
redshift as the H 2 component (see Fig. [8). A velocity shift be- 
tween H 2 and Ci of only 1.1 kms 1 is seen (see Table |3). 
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Fig. 2. Portion of the normalized spectrum of Q 2343+125. The resolving power is R 45 000. The best fit to the H2 lines is 
superimposed on the spectrum. Transition names are written AxBz, with A = L or W for the Lyman or Werner band systems; x is 
the quantum number of the vibrational upper state, while B is the branch designation (P, Q or R for AJ = +1,0, -1 respectively) 
and z is the rotational quantum number of the lower state. The fit parameters are given in Table [T] Zooms into individual lines are 
shown in small boxes that are shifted horizontally for clarity. Each box is 0.5 A wide. 



This supports the idea that neutral carbon and molecular hydro- 
gen absorptions originate from the same cloud. Thanks to the 
C 1 fine-structure level ratios, we can estimate the particle den- 
sity in the correspon ding cloud. Using a simple popratio model 
dSilva & V iegas 2001), taking into account the cosmic back- 
ground radiation at redshift z ~ 2 and using kinetic temperatures 
between T = 1 88 K and T = 300 K and a radiation field intensity 
between 1 and 10 times the Galactic value, leads to densities be- 
tween 20 and 40 ctrT 3 . A higher UV flux would result in a lower 
density, whereas a lower temperature would lead to a higher den- 
sity, that would be about twice higher t han derived above fo r the 
mean Galactic temperature (T ~ 77 K; lRachford et al.ll2002l) . 

4.3. Metal content 

Ten highly blended components are needed to reproduce the 
whole absorption profile of low ionization species that is spread 



over about 200 km s _1 (see Fig. [9]). The profile is characterized 
by a central and strong component which is the only one detected 
in H2 and C 1. S 11 and Si 11 are detected along the whole profile. 
The decomposition is well constrained by the Fe 11 A 1608 profile 
that is redshifted in a region with high SNR and is nowhere sat- 
urated except for the central component which is constrained by 
the Fe n/l/12249,2260 transition lines. Crn, Znn, Nin and Mgi 
are detected in, respectively, seven, eight, eight and three com- 
ponents. Ni is also detected over the whole profile, but because 
of many blends with intervening Ly-a absorptions, we consider 
the corresponding column densities as upper limits only, except 
for the central component where the column density is derived 
thanks to the three detected N 1 transitions. Mg 11 is also detected 
in the central component and possibly also in component # 9 (see 
Table 2]i. We nevertheless consider the Mgn column density in 
component # 9 as an upper limit because of a possible blend with 
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Table 2. Low-ionization metals column densities in the z a bs = 
2.431 DLA toward Q 2343 +125. 
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Fig. 5. Absorption profiles of low-ionization and neutral species 
in the damped Lyman-o- system at z a b s = 2.431 toward 
Q 2343+125. The dotted profiles in some panels mark absorp- 
tion from another element. 



Ly-a forest absorptions. Upper limits at 3 cr are derived for non- 
detected components. 

Some additional absorption is present in the blue and red 
wings of the Fen /U608 profile. Unfortunately there is no 
stronger Fen absorption line in our spectrum to constrain the 
column densities well. If we include these wings, the total Fe n 
column density is increased by about 7%. 

The result of Voigt-profile fitting is shown on Fig. [9] and 
summarized in Table [4] Summing up the column densities mea- 
sured in all components, we derive a metallicity of [Zn/H] = 
-0.85 ± 0.06. 



# 


'abs 


Ion (X) 


logiV(X) 


b 


Av-x/H, 










[km s- 1 ] 


[kms- 1 ] 


1 


2.43105(7) 


Ni 


13.49 + 0.02 


9.3 ±1.0 


_ 






O) 


15.46 ± 0.02° 










Mgi 


< 11.55* 










Sin 


14.31 ±0.08 










Pn 


12.16 ±0.09 










Sn 


13.41 ±0.13 










An 


11.64 ±0.26 










Cm 


12.00 ± 0.09 










Mnii 


11.51 ±0.10 










Fen 


13.57 ±0.08 










Ni [i 


12.39 ±0.14 










Zn ii 


< 11.00* 






2 


2.43116(9) 


Ni 


14.20 ± 0.01 


5.1 ±0.4 








Oi 


> 15.72 c 










Mg i 


11.94 ±0.06 










Sin 


14.55 ± 0.07 










Pn 


12.53 ±0.06 










Sn 


14.20 ± 0.04 










An 


12.82 ±0.02 










Cm 


12.31 ±0.06 










Mnii 


11.67 ±0.08 










Fe ii 


13.96 ±0.05 










Ni ii 


12.84 ±0.06 










Zn ii 


11.71 ±0.04 






3 


2.43129(3) 


Ni 


14.19 ±0.01 


5.6 ±0.5 


+0.8 






Oi 


> 15.92 c 










Mg i 


11.87 ±0.08 










Si n 


14.70 ± 0.04 










Pn 


12.64 ± 0.04 










Sn 


14.24 ± 0.04 










An 


12.78 ± 0.02 










Cm 


12.34 ± 0.05 










Mnii 


11.89 ±0.04 










Fe ii 


14.00 ± 0.04 










Nin 


12.97 ± 0.04 










Znn 


11.74 ±0.04 
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2.43143(8) 


Ni 


13.78 ±0.01 


5.6 ±0.5 








Oi 


> 15.81 c 










Mg i 


1 1.57 ± 0.13 










Sin 


14.45 ± 0.04 










Pn 


12.22 ± 0.06 










Sn 


13.88 ±0.04 
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12.32 ±0.04 










Cm 


12.14 ±0.05 










Mn ii 


11.66 ±0.05 










Fen 


13.90 ±0.03 










Nin 


12.71 ±0.05 










Znn 


11.57 ±0.04 
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2.43158(1) 


Ni 


13.18 ±0.03 


5.8 ±0.4 


— 






Oi 


14.82 ± 0.06° 










Mgi 


< 11.55° 










Sin 


13.80 ±0.09 










Pu 


< 11.85° 










Sn 


13.33 ±0.08 
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11.69 ±0.17 










Cm 


11.92 ±0.05 










Mnii 


11.28 ±0.09 










Fe ii 


13.40 ±0.04 










Nil! 


12.45 ± 0.06 










Zn ii 


< 11.00° 







" O i has been fitted using fixed Doppler parameters previously derived 
from the other ions. Errors on yV(0 1) could therefore be underestimated. 
b Upper-limits are given for non-detections at the 3 cr confidence level. 
c N(0 1) is considered as a lower-limit in the three saturated components 
because of uncertainties in the background subtraction. 



5. The system at z abs = 2.402 toward HE 0027-1 836 

The very broad Ly-a absorption line at z aDS = 2.402 lies in 
the wing of the Ly-a emission line from the quasar (z em = 
2.56) where the normalization of the spectrum is problematic. 
We therefore fitted a Voigt-profile to the Hi Ly-/3 line (see 



Fig. [TOl to derive a total neutral hydrogen column density of 
log N(Hi) = 21.75+0.10. This is one of the highest column den - 
sities measured in QSO-DLAs (see, e.g. lProchaska et ai1l2005l) . 
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Table 4. Low-ionization metals column densities in the z a b s = 1.989 DLA toward HE 23 18- 1107. 
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Znn 


< 10.80 








Znn 


10.99 ±0.13 






2 


1.98827(8) 


Ni 


< 13.79 


8.6 ± 0.2 


— 7 


1.98918(3) 


Ni 


< 14.09 


9.5 ± 1.6 


— 






Mgi 


< 11.50 




— 




Mgi 


< 11.50 




— 






Mgn 


blend 








Mgn 


< 14.00 










Sin 


14.05 ± 0.02 








Sin 


14.31 ±0.02 










Sn 


13.75 ± 0.02 








Sn 


14.18 ±0.02 










Cm 


11.96 ±0.08 








Cm 


11.83 ±0.11 










Fen 


13.71 ± 0.04 








Fen 


13.92 ± 0.08 










Nin 


12.58 ± 0.05 








Nin 


12.90 ± 0.02 










Znn 


11.18+0.12 








Znn 


11.63 ±0.04 






3 


1.98845(3) 


Ni 


< 13.53 


8.4+1.5 


— 8 


1.98942(1) 


Ni 


< 14.01 


10.5 ± 1.8 


— 






Mgi 


< 11.50 




— 




Mgi 


11.80 ±0.13 




— 






Mgn 


blend 








Mgn 


< 14.00 










Sin 


14.23 ± 0.05 








Sin 


14.40 ± 0.03 










Sn 


14.00 ± 0.04 








Sn 


14.07 ± 0.06 










Cm 


< 11.70 








Cm 


12.17 ±0.05 










Fen 


13.43 ± 0.03 








Fen 


14.03 ± 0.03 










Nin 


12.63 ± 0.04 








Nin 


12.96 ± 0.02 










Znn 


11.53 + 0.05 








Znn 


11.44 ±0.13 






4 


1.98865(0) 


Ni 


< 13.62 


10.7 ± 4.0 


— 9 


1.98956(4) 


Ni 


< 13.39 


1.6 ± 1.0 


— 






Mgi 


< 11.50 




— 




Mgi 


< 11.50 




— 






Mgn 


blend 








Mgn 


< 14.48 










Sin 


14.62 ± 0.05 








Sin 


13.85 ± 0.08 










Sn 


14.43 ± 0.05 








Sn 


13.61 ±0.08 










Cm 


12.37 ± 0.03 








Cm 


< 11.70 










Fen 


14.13 ±0.02 








Fen 


13.42 ± 0.04 










Nin 


13.11 ±0.02 








Nin 


< 12.30 










Znn 


11.72 ±0.04 








Znn 


< 10.80 






5 


1.98888(9) 


Ni 


14.55 ± 0.02 


4.2 ± 0.3 


+0.7 10 


1.98970(7) 


Ni 


< 14.01 


9.2 ± 1.9 


— 






Mgi 


12.32 ± 0.05 




-0.4 




Mgi 


12.04 ± 0.08 










Mgn 


14.94 ± 0.04 








Mgn 


blend 










Sin 


14.72 ± 0.01 








Sin 


14.42 ± 0.03 










Sn 


14.50 ± 0.04 








Sn 


13.95 ± 0.04 










Cm 


12.60 ± 0.01 








Cm 


12.26 ± 0.04 










Fen 


14.35 ± 0.01 








Fen 


13.99 ± 0.08 










Nin 


13.18 ±0.01 








Nin 


12.84 ± 0.03 










Znn 


11.97 ±0.01 








Znn 


11.53 ±0.05 







" 3 cr upper-limits are derived for non detected components, and marked by "<". 

* Components for which the column density is uncertain because of possible blends are considered as upper-limits, and marked by "<" 



5.1. Molecular content 

The wavelength range of the UVES spectrum allows us to ob- 
serve a large number of H2 absorption lines at the redshift of 
the DLA. A single H2 component is detected in the rotational 
levels J = to J = 5 with a possible detection of the J = 6 rota- 
tional level. In addition to the overall normalization of the spec- 
trum, the continuum has been carefully normalized over about 
100 to 800 km s locally around each line. It must be noted that 
this system is exceptional amongst H2-bearing DLAs, with more 
than 70 detected H2 transitions. The absorption can be fitted with 
a single component and most of the lines are not strongly satu- 
rated. In addition, at this redshift (z a bs ~ 2.4) blending is not 
strong as the Ly-a forest is not dense and the data are of high 



signal-to-noise ratio and high spectral resolution (R 50000). 
Therefore we are in an ideal position to derive with high con- 
fidence the column density N and Doppler parameter b(J) for 
each rotational level J. The fit was performed by x 2 minimiza- 
tion, considering different b- values for each rotational level. The 
relative optical depths of the different absorption lines with dif- 
ferent oscillator strengths provide strong constraints on N but 
more importantly on b. This is particularly true for rotational 
levels J = 2 and J = 3 for which a large number of transitions 
with different oscillator strengths are observed. Some of the ab- 
sorptions are in the intermediate regime where the equivalent 
widths strongly depend on b (corresponding to the flat part of 
the curve-of-growth). Note that damping wings appear for some 
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Fig. 6. Fit to the damped Lyman-*? line at z a bs = 1 .989 toward 
HE 23 18-1 107. The derived neutral hydrogen column density is 
log N{Hi) = 20.68 + 0.05. The dashed curves correspond to the 
best-fitted Voigt profile (black) and the associated errors (red). 
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Fig. 8. C i ground-state fine-structure lines from the z a bs = 1 .989 
DL A toward HE 23 1 8 - 1 1 07 . 
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Fig. 7. H2 absorption profile at z a b s = 1.98888 toward 
HE 2318-1 107. The model fit is over-plotted. 



Table 3. H2 and C 1 column densities in the z a b s = 1 .989 DLA 
toward HE 2318-1107. 



the J = 2 and J = 3 features. Two consistent absorption features 
are detected for rotational level J = 6, but the lines are too weak 
to perform a fit with an independent b. We therefore have chosen 
to fix b(J - 6) equal to b(J = 5) previously measured indepen- 
dently. The derived J = 6 column density should be preferably 
considered as an upper limit due to possible blends (see Table[5]l. 
The results of the Voigt-profile fitting are shown in Figs. [19] to 
|25l and summarized in Table [5] The molecular fraction in this 
system is log/ = -4.15 + 0.17. The important finding here is 
that the Doppler parameter b increases with J. 



Table 5. H2 and C 1 column densities in the z a b s = 2.402 DLA 
toward HE 0027-1836. 



r Zabs 


level 


logiV(H 2 , J) 6 


b 

[km s" 1 ] 


T -, c 
[K] 


3 2.40183(4) 


J = 


16.75 ± 0.06 


1.06 ±0.13 






J= 1 


17.15 ±0.07 


1.46 ± 0.14 


134 + 3' 7 




J = 2 


14.91 ± 0.02 


2.67 ± 0.08 


88 + » 




J = 3 


14.91 ± 0.01 


3.77 ± 0.07 


141 4 9 ° 




J = 4 


14.22 ± 0.01 


4.61 ± 0.25 


213 + !* 




J = 5 


14.02 ± 0.03 


6.17 ±0.69 


262 +lS 

— 16 




J = 6 


< 13.53 


6.17 rf 


< 375 



r 


Zabs 


level 


log A/(H 2 ,J)* 


b 

[kms" 1 ] 


[K] 


5 


1.98888(2) 


J = 
J= 1 
J = 2 
J = 3 


14.72 ± 0.03 
15.28 ±0.03 
14.84 ± 0.03 
< 14.7 rf 


3.6 ±0.2 


188 + ^ 8 
384^ 3 
<355 




r 


Zabs 


logMCi) 


logiV(Cf) 


b 

[km s" 1 ] 


Avci/H 2 

[km s" 1 ] 


5 


1.98889(3) 


12.63 ± 0.02 


12.30 ± 0.04 


2± 1 


+ 1.1 



8 The numbering refers to that of the metal components (see Table|4j. 
* The errors on the H2 column densities are the 1 a uncertainty from 
fitting the Voigt-profiles. 

c The excitation temperatures are calculated using the 3 <r range for H 2 . 
3 a upper-limit. 



lines of the J = and J = 1 rotational levels. For the highest ro- 
tational levels (J = 4 and J = 5), absorption lines are optically 
thin, making b slightly less constrained than for lower rotational 
levels. However, it is apparent that these lines are broader than 



r zabs logW(Ci) logAT(Cf) b Av Ci/ h 2 

[kms~'] [kms -1 ] 

3 2.40185(8) 12.25^-^ < 12.27 0.6-3.0 +2T~ 

" The numbering refers to that of the metal components (see Table[6l(. 
b The error on the H 2 column densities are the 1 cr uncertainty from fit- 
ting the Voigt-profiles. 

c The excitation temperatures are calculated using the 3cr range for 
log2V(H 2 ,D. 

d The Doppler parameter for J = 6 has been fixed to be equal to b(J = 5) 
derived independently. 



5.2. Neutral carbon 

Neutral carbon is detected in one component at about the same 
redshift as the H2 component. The absorption line is optically 
thin, and the column density therefore does not depend much 
on b. The corresponding C 1* absorption is weak and detected 
just above the 3cr confidence level (logA^(Cf) = 12.2). The 
derived column density is thus uncertain and should be consid- 
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Fig. 9. Absorption profiles of low-ionization and neutral species 
at z abs = 1.989 toward HE 23 18- 1107. Note that Mgi ,12026 
is partly blended with Znn /12026, and so is Znn ,12062 with 
Cr n /12062. N i lines are blended altogether as well as with Ly-a 
forest absorptions but the column density can be derived in the 
central component. 



ered as an upper limit: logA^(Cf) < 12.27 (see Table |5j. We 
observe a small velocity shift between the positions of the C i 
and H2 components but no velocity shift between C 1 and com- 
ponent #3 of the low-ionization metal profile. The upper limit 
on N(C f)/N(C 1) leads to an upper limit on the neutral hy- 
drogen density, «h < 150 cirT 3 , for a typical temperature of 
100 < T < 200 K and a Galactic UV radiation field (see also 
iNoterdaeme et al.ll2007l) . 




-aooo 3000 

Relative velocity (km s _1 ) 

Fig. 10. Fit to the Ly-/3 line at z abs = 2.402 toward 
HE 0027- 1836. The neutral hydrogen column density is 
logTY(Hi) = 21.75 + 0.10. The dashed curves correspond to the 
best-fitted Voigt profile (black) and the associated errors (red). 
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Fig. 11. Neutral carbon absorption lines toward HE 0027-1836. 
The C i* absorption is weak and the derived column density is 
thus considered as an upper limit. 



z 5.3. Metal content 



The profile is made of two strong well separated features (Av ~ 
25 kms -1 ), the first one being a blend of two components. The 
second feature, corresponding to the third and also the strongest 
component in the profile , is associated with the H2 component 
but is significantly broader. We detect and measure column den- 
sities for a large number of species, namely N 1, Mg 1, Mg 11, Si 11, 
P 11, S 11, Ar 1, Ti 11, Cr 11, Mn 11, Fe n, Ni n, and Zn n. In addition, 
we measure 3 cr upper limits on the column densities of Pb 11 
(resp. Cu 11) for each component: log A^(Pb 11) < 11.1 (respec- 
tively, log N(Cu 11) < 1 2. 1 ) from the non-detection of Pb 11 A 143 3 
at SNR ~ 46 (resp. Cun ,11358; SNR ~ 16). 

To derive the characteristics of the lines, we first determined 
the ^-values by fitting unblended and good SNR lines. We then 
fixed the obtained b values and added all other detected tran- 
sitions. Although relative wavelengths were relaxed when fit- 
ting the lines, the final fit shows no velocity-shift between the 
Ci component and the third metal component. Note that the 
best oscillator strength found in the literature for Fen /ll 112, 
/ = 0.00629 ( jHowket alj|2000h . seems to be underestimated 
(see Fig. IT2T> unless there is extra absorption (possibly from the 
Ly-q forest) a t that wavelength. The oscillator strength from 
iMortonl d2003l) is even smaller, / = 0.00446. Usin g the mea- 
sured i vac = 1253.805 for Snil253 dMortonll2003l) induces a 
small but significant shift between the observed and the synthetic 



profiles. To make th e two spectra consistent the Ritz wavelength 
A vac = 1253.811 dMortonlll99lb had to be used instead. The 
metallicity derived for this system is [Zn/H] = -1.63 + 0.10. 
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Table 6. Low-ionization metals column densities in the z a b s = 
2.402 DLA toward HE 0027-1836. 



Table 7. C n* column densities in the z abs = 2.402 DLA toward 
HE 0027-1836. 



# 




Ion (X) 


logJV(X) 


b 

[kms" 1 ] 


Av X /h 2 
Avx/ci 
[km s~>] 


1 


2.40150(4) 


Ni 


14.53 ± 0.24 


8.4 ± 1.5 


— 






Mgi 


11.52 ±0.09 




— 






Mgn 


15.32 ± 0.14 










Sin 


14.84 ±0.11 










Pn 


12.62 ± 0.30 










Sn 


14.41 ± 0.10 










Ari 


13.64 ± 0.06 










Tin 


11.87 ±0.29 










Cm 


12.71 ± 0.04 










Mnn 


11.93 ±0.16 


" 








Fen 


14.36 ± 0.05 










Nin 


12.88 ± 0.10 










Znn 


11.99 ±0.10 






2 


2.40159(6) 


Ni 


14.51 ±0.22 


4.7 ± 0.9 


— 






Mgi 


11.96 ±0.09 




— 






Mgn 


15.14 ± 0.17 










Sin 


15.06 ± 0.07 










Pn 


12.50 ± 0.30 


" 








Sn 


14.68 ± 0.03 










An 


13.64 ± 0.06 










Tin 


11.90 ±0.11 










Cm 


12.81 ± 0.04 


" 








Mnn 


12.31 ±0.08 










Fen 


14.39 ± 0.05 










Nin 


13.21 ± 0.05 










Znn 


12.09 ± 0.06 


" 




3 


2.40185(8) 


Ni 


15.04 ± 0.09 


4.7 ± 0.2 


+2.1 






Mgi 


12.41 ± 0.05 




0.0 






Mg ii 


15.81 ±0.02 










Sin 


15.45 ± 0.02 










Pn 


12.87 ± 0.40 










Sn 


14.98 ± 0.03 










An 


14.24 ± 0.02 










Tin 


12.41 ± 0.03 










Cm 


13.08 ± 0.01 










Mnn 


12.60 ± 0.02 










Fen 


14.66 ± 0.03 










Nin 


13.42 ± 0.02 










Znn 


12.60 ± 0.01 







A study of the highly ionized species is presented by lFox et"aTI 



A study o 
(2007a,b) 



5.4. Cn* and the ambient UV flux 

Strong and saturated C n*/H335 absorption is seen in this sys- 
tem (see Fig. [T3l ). We can use however the weaker transition 
at Aresi — 1037A to derive reliable column density estimates 
(see Table |7). During the fit redshifts of the C n* components 
were constrained to be the same as those of other low-ionization 
species. In turn, Doppler parameters could not be given the same 
values. In any case the ^-values for C n* are not well constrained, 
especially for the weak component # 1 but this does not strongly 
affect the derived column densities. The associated uncertain- 
ties were derived by considering a range of b values. A large 
Doppler parameter is needed for component #3. This probably 
reveals the multipl e component nature of the feature. Following 
IWolfe et alJ (|2003), we can estimate the UV ambient flux in 
the system, considering the equilibrium between photo-electric 



r 




log/Y(Cn*) 


b 








[kms" 1 ] 


l 


2.40150(4) 


12.65 ± 0.30 


3-8 


2 


2.40159(6) 


13.80 ±0.10 


4-7 


3 


2.40185(8) 


14.15 ±0.25 


7-9 



" The numbering refers to that of the low-ionization metal components 
(see Table |6). 



heating and cooling by C n 158 yum emission. The heating of the 
gas by photo-electric effect on dust grains is proportional to the 
dust-to-gas ratio, (* = 10 [Zn/H] (l - 10F e/Zn] )), and to the UV flux, 
F\jy. The cooling rate is pro portional to the C n 15 8 //m emission 
per H atom, as estimated bv lPottasch et alJ (Q979): 



N(Cu*)hv u ,A ul _j 

L = erg s per H atom, 

iVTHl) 5 F 



(1) 



where A u \ is the coefficient for spontaneous emission of the 
2 ^3/2 — * 2 Pi/2 transition and hv u i the corresponding energy. We 
find \ogl c — -27 ± 0.2. The conditions in the gas can be com- 
pared to that in the Galactic ISM, considering the same type of 
dust grains for both media: 



1 



I* 



(2) 



From Fig. 2 of IWolfe et all d2003l) (see also lLehner et~ai1l2004h 
it can be seen that log if 1 « -26.6 in the Galaxy at N(H i) = 
10 21 - 75 cm" 2 . We then estimate F uv IF g ^' v ~ 20 for k a 0.02. The 
UV flux in the DLA is about an order of magnitude higher than 
in the Galactic disk. Note that we ignored the dependence of 
the photo-electric heating on the temperature and density of the 
gas (jW em gartner & Drainel200ll) . and that the result is therefore 
only valid for temperatures and densities similar to that in our 
Galaxy. 




100 



50 100 
Relative velocity (km s -1 ) 

Fig. 13. C n* absorption profile in the DLA system toward 
HE 0027- 1836. The absorption seen at v - +35 kms 1 in the 
first panel is due to the H2L5-ORI transition. 



6. Excitation of H 2 toward HE 0027-1 836 

In this Section, we analyze in detail the excitation of H2 at 
z a b s = 2.40183 toward HE 0027-1836 where absorptions from 
rotational levels J = 1 to J = 5 and possibly J = 6 are seen in a 
single component. 
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Fig. 12. Absorption profiles of low-ionization and neutral species from the z a b s = 2.402 DLA system toward HE 0027-1836. The 
synthetic best-fit spectrum is superimposed to the observed spectrum. The dashed lines in some panels mark absorption from another 
element. For example, there is some MgiT2026 absorption near the ZnnT2026 one. The transitions Feii/i/i/il063, 1064,1112 are 
shown but were not used during the fitting process because of uncertainties in atomic parameters. 



6.1. Excitation temperatures 

Figure [14] shows the H2 excitation diagram in this system. The 
graph gives for the different rotational levels 7V(H2,J)/gj, where 
gj is the statistical weight of the level, versus the relative energy 
between that level and J = 0. The slope s of a straight line in the 
graph is inversely proportional to the excitation temperature, s = 
-1/ (Xi-y In 10). The kinetic temperature is generally estimated 
by 7Vi assuming the J = 1 level is thermalised. Here, 7Vi - 
134 K. 

The H2 excitation diagram is characterized by two excitation 
temperatures, 7Vj ~ 90 K for J < 2 and T 2 -j ~ 500 K for J > 2 
(see Fig.[T4"l>. Therefore, other processes than collisions at a tem- 
perature of the order of T ~ 100 K (the approximate kinetic tem- 
perature) are at play. Formation pumping (just after formation, 
H2 molecules are in high-J levels) and UV pumping from low-J 
levels upwards are likely to be the most important of these pro- 
cesses. Assuming these are the main processes allows us to esti- 
mate the strength of the surrounding UV radiation field. Writing 



the equilibrium between the processes popul ating and depopu- 
lating the J = 4 rotational level [see Eq . 5 in [No terdaeme et alj 
d2007l) . also iHirashita & Ferraral d2005b : ICui et alj d2005h l. we 
can estimate the photo-absorption rate in the J = level, fio, in- 
side the H2-bearing cloud:j6o 2.5xl0~ n s _1 . This is a very low 
value, probably due to self- and dust-shielding by the outer lay- 
ers ofJl^^l^u^TT^corrective term for shielding is expressed 
as dDraine & Bertoldilll996l) : 



S shi 



shield 



N(H 2 ) 
10 14 cm- 2 



-0.75 



(3) 



This is the product of two terms, the first term, 5h,, is due to self- 
shielding, the second, S dust, is due to dust extinction. The total 
H2 column density in the cloud is A^(H2) = 10 17 3 cm 2 , leading 
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Fig. 14. Excitation diagram of H 2 at z abs = 2.40183 
HE 0027-1836. It is apparent that it is characterized by 
citation temperatures: Tq-j ~ 90 K for J < 2 and 7Vj ~ 
forJ>2. 



4000 

toward 
two ex- 
- 500 K 



to Sh-. — 0.00335. Dust extinction is calculated using the dust 
optical depth, tuv: 



tuv = 0.879 



0.1pm/ \2gcm 



(-U( N » ) (4) 

\10- 2 /\10 21 cm- 2 / V ' 



where a is the radius of a grain, g i s the dust-to-gas mass rati o 
and 6 is the grain material density ( Hiras hita & Ferra ra 2005). 
This optical depth is difficult to estimate, due to the unknown 
type of dust and the uncertainty on the exact H i column density 
in the H2-bearing component. We can assume a = 0.1 /urn, 6 
2 g ctrT 3 , and scale the dust-to-gas mass ratio g with the dust-to- 
gas ratio k = 10 [Zn/H] (l - 10 [Fe/Znl ) dProchaska & Wolfell2002l) . 
We then get: 



or equivalently, 



r uv = 0.879 



N(Zn) JV(Fe) 



10[ Zn / H ]o iq[^/h] / \ 10 21 cm 



1 



(5) 



(6) 



where the column densities are those measured in the H2 -bearing 
component. We obtain S dust = 0.94 and S shield = S h 2 S dust - 
0.003. Note that the self-shielding dominates the total shield- 
ing. We can then estim ate the UV flux outside the cloud (see 
iNoterdaeme et alj|2007l) : 

y LW * 8xl0 '" A) * 6.3 x lO^ergs- 1 cm- 2 Hz" 1 sr" 1 . (7) 

S shield 

This is about 20 times higher than in the solar v icinity (7lw,q — 
3.2 x 10 _20 ergs _1 cm" 2 Hz" 1 sr ': lHabing1ll968h . Again, all this 
assumes that the high rotational levels of H2 are populated by 
formation- and UV-pumping. If collisional excitation, in case of 
turbulence or shocks for example, is playing an important role, 
then the derived UV flux should be considered as an upper limit. 

6.2. Velocity dispersion 

In Fig. Q3] we plot the Doppler parameter of the lines from a 
given rotational level, b, as a function of the energy of the ro- 
tational level J. It is apparent that higher J levels have broader 



lines. This effect has already been observed in the local inter- 
stellar medium in some cases. It has first been derived from 
curve-of-growth analysis in several Coper nicus observations 
dSpitzer & Cochranlll973HSpitzer et alJI 19741) . then observed di- 
rectlv bv iJenkins & Peimbertl (ll997T with R = 120000 data 
from the Interstellar Medium Absorption Profile Spectrograph 
(IMAPS). In that case however the increase in the line broaden- 
ing is associated with a regular shift i n velocity. More re cently, 
the same effect has been reported bv lLacour et al.l d2005l) along 
four lines o f sight toward early- type Galactic stars. In high red- 
shift DLAs. lLedoux et all (120031) already suggested that the lines 
of J > 2 rotational levels require higher b values than the lower 
J-level lines. It is however the first time that a systematic effect 
is seen beyond any doubt at high redshift. 
In case the broadening is only thermal, the Doppler parameter is 
related to the kinetic energy by: 



b(T) = y/2k B E k (J)/(3m H ) 



(8) 



If we assume that the excess of kinetic energy of molecules in 
some J level compared to the kinetic energy of molecules in the 
J = level is directly proportional to the energy of the rotational 
level, i.e. ^(J) - E^Q) = a Eq-j, with a a constant, then we can 
try to fit the observed trend with the expression: 



b(J) = V2M£*(0) + a£o-j)/(3m H ) 



(9) 



where ^(O) is the kinetic energy of rotational level J = 0, for 
which we use the excitation temperature 7o-i - 134 K previ- 
ously determined. The least square minimization gives Ef,(J) - 
Ek{0) + 2.34 £o-j. The result is shown on Fig. [15] One can see 
that this very simple assumption fits the data very well. It is not 
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Fig. 15. The Doppler parameter b as a function of the energy 
between the rotational levels (from J = to J = 5) in the DLA 
system toward HE 0027-1836. The smallest errors are found for 
J = 2 and 3, for which numerous lines with good SNR and in 
different optical depth regimes are observed. The dotted curve 
corresponds to a fit with a single free parameter (Eq.|9]l. 



easy t o find an explanation to this effect. IJenkins & Peimbertl 
dl997h argued that such trend in b seen along the line of sight 
to ^fOriA could be explained if the gas is located behind a J- 
shock. However, the column densities along the present line of 
sight are much larger and, more importantly, this would produce 
coherent velocity shifts between the lines from different J levels, 
which we do not observe. 
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UV pumping cannot produce different b values because 
the cascade following the UV absorption releases the energy 
through infrared photons, and such process does not result in 
any change in the molecule's kinetic energy. 

lLacour et al. I (120051) show that such behavior can be ex- 
plained by the consequences of ^-formation onto dust grains. 
After formation, the molecules escape with an excess kinetic 
energy left-over from the formation process. The highest rota- 
tional levels have a radiative life-time considerably shorter than 
the lowest J levels. This implies that low-J molecules have more 
time than hig h-J molecules to cool down through collision s with 
Hi (see, e.g. JSpitzer & Cochranll 19731; lLacour et alj|2005l) . This 
implies that, if other processes are negligible, the kinetic energy 
of high-J molecules coul d be larger th a n th at of low-J ones. It 
can be seen on Fig. 6 of lLacour et al.l (120051) that their model 
approximately reproduces the values measured here. However, 
their model requires a formatio n rate ~10 times higher tha n that 
measured in the Galactic ISM dJurall 1 975b iGrv et al.1120021 It is 
clear that at the low metallicity measured in the present DLA 
system, [Zn/H] = -1.63, this assumption can be rejected. 

Another explanation could be that the cloud is composed of 
several layers with a gradient of temperature. External layers 
would be warm and exposed to strong external UV flux. In these 
layers, the excitation of H2 could be large and mostly due to UV 
pumping. They would contribute mostly to the column densities 
of high J-levels. The internal layers would be cold and shielded 
from the external UV flux. It must be realized that the increased 
b with higher J is coupled with an excitation diagram described 
by two excitation temperatures, 7Vj ~ 90 K for J = 0-2 and 
T2-} ~ 500 K for J = 3-6 (see Fig.[T4"li. This kind of diagram has 
been explained in the Galactic ISM by invoking the association 
of a diffus e cold cloud with a warm Photo-Dissociation Region 
(PDR), see lBoisse et al.l(l2005l) . 

6.3. Models 

We used the Meudon PDR code dLe Petit et alj l2006t 
iGoicoechea & Le Bourlotl 120071) to model the system toward 
HE 0027- 1836 where H2 is detected in rotational levels J = to 
J = 5. This PDR model assumes a stationary plan-parallel slab 
of dust and gas of constant hydrogen density «h, illuminated by 
a UV radiation field and solves the radiative transfer, chemistry 
and thermal balance. The code used in this paper is a slightly 
modified version compared to the online onqj . We adopt th e 
Interstellar Radiation Field (ISRF) as given by iDraina (Q978), 
scaled by a factory. Note that the energy between 912 and 2400 
A for the Draine ISRF is 1.78 times that of the Habing radia- 
tion field. The grai n size distribution is assumed to follow the 
Math is et al.l d 1977b law with radii between 0. 1 and 0.3 fim. The 
mean Galactic dust extinction curve is used. Elemental abun- 
dances for C, N, O and S are scaled from solar abundances, 
[X/H] Q , of lMortonl (2003), using the observed sulfur metallic- 
ity, [S/H] (see Table[6]). Savage & Sembachl dl996l) have shown, 
in the context of diffuse Galactic clouds, that depletion varies 
significantly from one line of sight to another. Here, we test two 
extreme assumptions to fix the gas phase abundances of C and 
O: the same depletion on dust a s in the cool medium towards 
£ Oph (Savage & Sembach 1996) or no depletion. This gives in 
the first case C/H = 2.35 x 10~ 6 and O/H = 3.9 x 10~ 6 and in 
the second case 6.2 x 10~ 6 and 1.0 x 10~ 5 , respectively. CMB 
temperature is assumed to be 9.2 K (z — 2.4). 
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We built two grids of models: one corresponding to isother- 
mal models with r gas = 130 K corresponding to the observed 
value of 7Vi and one corresponding to models solving the ther- 
mal balance equations. For each grid we compute models with 
X = 1, 10, 20, 30 and n H = 10, 100, 150, 200, 400, 500 ctrT 3 . 
We checked the influence of the dust-to-gas mass ratio, g, or 
equivalently Cd = Nn/E(B-V), the ratio of total neutral hydro- 
gen column density to color excess, using g _ 0.01 and Cp = 
5.8 x 10 21 (corresponding to the ISM values; lBohlin et al.ll978h . 
0.001, 5.8 x 10 22 and 0.0001, 5.8 x 10 23 cm" 2 mag" 1 . Column 
densities are computed perpendicular to the plane parallel slab. 
Since the size of the cloud is not known, all models assume 
a semi-infinite geometry up to reproduce half of the observed 
NQI2). Then column densities are multiplied by 2 to simulate a 
slab of gas illuminated on both sides. The model is constrained 
by the condition that the total observed NQI2) should be repro- 
duced. 

The best iso-thermal models are obtained with g = 0.001 and 
Co = 5.8 x 10 22 cm^mag -1 corresponding, as expected, to a 
low dust content (note that if we scale depletion with the ob- 
served [Zn/Cr] ratio, see previous Section, we obtain g =0.002). 
Figure[T6lpresents the excitation diagrams obtained with isother- 
mal assumption in two cases. The best fit to the data for all lev- 
els is obtained with «h = 150 crrT 3 , % — 10 (about 20 times the 
Habing radiation field). This model reproduces reasonably well 
all column densities except that of the J = 2 level which is over- 
estimated by nearly an order of magnitude. The low observed 
value of N(i-2) relative to other column densities is the special 
feature of this DLA system. In this model, the molecular fraction 
is / = 10~ 4 so H is the main collisional partner for H2. It must 
be noted that H/H2 collision rates are poorly known mainly be- 
cause of reactive collisi ons. However using c ollision rates from 
iFlower & Rouefj dl998l) instead of those from lMandv & M artin 
d 19931) does not change the results significantly. Amongst mod- 
els of total H2 column density corresponding to the observed 
one, only those with low-density, «h - lOctrT 3 , and weak radi- 
ation field, x — 1 (Fig- fl6t , can reproduce the observed iV(H2, 
J=2). Higher densities tend to thermalize J = 2 whereas higher 
radiation fields tend to pump this level too much. The influence 
of the distribution of H2 amongst ro-vibrational levels after for- 
mation on dust has also been tested. It can be noted that a few 
lines of sight with 7o_2 < 7Vi have been observed by FUSE 
and Copernicus in our Galaxy or the Magel lanic clouds. Th e 
HD 108927 line-of-sight has been modelled bv lGrv et alJd2002l) . 
The authors reached similar conclusions. However, with such 
low density and UV flux («h — 10 ctrT 3 , x — 1)> populations 
of J > 2 levels will be under-estimated by more than one order 
of magnitude. Note, however, that in these models, N{C 1) and 
N(Cf) are in good agreement with observations. Depending on 
the adopted depletion of C into dust grains, N(C 11*), N(C 1) and 
N(C 1*) are reproduced within a factor of 2 to 3 which is accept- 
able considering the uncertainties in the model. N(S 11), iV(Mg 1) 
and N(Mg 11) are also in good agreement with observations. 

Thermal balance models lead to similar conclusions. These 
models are, in principle, more constrained as they have to re- 
produce 7o-i. However, they are dependent on assumptions on 
(i) the dust composition since the photo-electric effect on dust is 
the main heating mechanism and (ii) the elemental abundance 
of oxygen which controls cooling. Best model is again for g 
= 0.001 and C D = 5.8 X 10 22 cm 2 mag 1 . As for isothermal 
models, a low density is required to avoid to overestimate iV(H2, 
J=2). With the first set of abundances (same depletion as in the 
cool medium towards f Oph), the best model is «h = 50 ctrT 3 
and x - 1- In this case, 7Vi = 90 K. A lower density of 10 
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Excitation diagrams obtained for isothermal models 
130 K. Grain parameters are g = 0.001 and Co = 



However, the probability that the line of sight crosses such 
a clump is small especially when the model requires that this 
small clump must be embedded in a particularly strong UV ra- 
diation field. Note also that there is no velocity shift between the 
absorption lines of the different J levels suggesting that they are 
not produced in very different locations. Another possibility to 
explain the e xcitation is turbulen t dissipation either as vortices 
or C-shocks dJoulain et all 1 19981: ICecchi-Pestellini et al.1 [2005 ; 
iGredel et all 120021: iLe Petit et al.ll2004l) . In that case, increased 
temperature in turbulent vortex would be responsi ble for the ex- 
citatio n of the higher H2 rotational levels. Indeed, Joulair Tet alj 
(1998) show that there should be no shift between the cen- 
troids of different species and that the widths of the lines should 
be larger for higher excitatio n as observed in the present case. 
ICecc hi-Pestelli ni et al.l (12005b show that a small amount of hot 
gas located in turbulent dissipative cells can explain the H2 exci- 
tation. However, temperature excitation as high as 7Vj = 500 K 
for a low total N(H2, J>2) may be difficult to be reproduced by 
such models. This assumption probably deserves more detailed 
investigation. 
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Fig. 17. Excitation diagrams obtained with thermal balance 
models. Grains parameters are g = 0.001 and Co = 5.8 x 
10 22 cirT 2 mag -1 . 



cm increases 7Vi to 180 K because of the low abundance of 
coolants. With the second set of abundances (no depletion), n H 
= 10 cirT 3 and ^ = 1 is favored and gives 7o_i = 96 K. The 
corresponding excitation diagrams are presented in Fig. [171 

In the two models, 7o-i is representative of the kinetic tem- 
perature which varies from about 100 to 88 K when the depth 
in the cloud increases. The best-fit model with a depletion fol- 
lowing that towards £ Oph gives N(Cu") = 5 x 10 13 , N(Ci) = 
8.0 x 10 11 and iV(Ci*) = 7.8 x 10 11 cm 2 . The one with no de- 
pletion gives 1.4 x 10 14 , 5.8 x 10 12 and 3.2 x 10 12 cirT 2 for the 
same species. Both models are within a factor of two from ob- 
servations (see Tables [7] and 0. 

The population of high-J levels cannot be reproduced by 
the above low-density models and a PDR component has to be 
added. High density and strong radiation field are required in 
order to explain H2 excitation by fluorescence. Fig. [T7]gives an 
example in which the high-J populations are produced in a clump 
of density «h = 500 cm -3 with x — 30. The corresponding size 
is 0.3 pc. Note that the value of the UV radiation field intensity 
is consistent with what is estimated analytically in Sects. l6"Tl and 



7. Comments on metallicity and depletion 

In Fig. [18] we compare the depletion patterns for the three sys- 
tems presented here to the typical depletions observed in cold 
and warm gas of the Galactic disk and gas in the Galactic 
halo. The different points correspond to the different compo- 
nents in each system in which Znn is detected. It is appar- 
ent that the Galactic halo depletion pattern represents best the 
observed abundances. The differences between components are 
small. In the DLA system toward HE 0027-1836, the depletion 
found in the ^-component is the highest in the system how- 
ever: [Zn/Cr] = 0.54, while it is ~ 0.3 in the other components. 
Similarly, [Zn/Fe] = 0.77 in the H2-bearing component, while 
[Zn/Fe] = 0.46 and 0.53 in the other two components. Note 
that silicon is almost non-depleted in all components of the three 
systems. This is similar to what is obser ved through some lines 
of sight in the Small Magellanic Cloud dWeltv et al.ll200ll) . The 
small depletion factor, together with low metallicities, implies 
that the dust content is small and can explain the low observed 
molecular fraction. The fact that we do not detect H2 in other 
components cannot only be due to the little difference in dust-to- 
gas ratio, resulting in a lower H2 formation rate. The main reason 
is probably the lower column densities in the components (espe- 
cially for HE 2318-1 107). We summarize the total abundances 
found in the three systems in Table [8] 

8. Conclusion 

We have presented a detailed analysis of three H2-bearing 
damped Lyman-o' systems. Two of them are reported here for 
the first time. This brings the number of known high-redshift 
(Zabs > 1-8) H2-bearing DLAs to twelve. All three systems 
have low-molecular fractions (log / < -4). Only one DLA sys- 
tem with such a low molecular fr action was reported be fore, at 
z abs = 3.025 toward Q 0347-383 dLevshakov et al.ll2002l) . 

The depletion patterns observed in the components of the 
three systems are very similar to what is observed in gas located 
in the Galactic halo, probably because of similarly low metal- 
licities. The depletion is not very different, although slightly 
larger, in the H2-bearing components compared to the other 
components in the systems. This is different from what is seen 
in systems with larger molecular fractions where large deple- 
tion factors are usually observed in H2-bearing components 
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Table 8. Summary of metals abundances in newly discovered H2-bearing DLAs. 



QSO 


Q 2343+ 125, Zem = 2.51 


HE 2318-1 107, z em = 2.96 


HE 0027- 1836, Zem = 


: 2.56 


Zabs 


2.431 


1.989 




2.402 




logyV(H i) 


20.40 


±0.07 


20.68 


±0.05 




21.75 


±0.10 




logjV(H 2 ) 


13.69 


±0.09 


15.49 


±0.03 




17.30 ± 0.07 




log/ 


-6.41 


±0.16 


-4.89 


± 0.08 




-4.15 


±0.17 




Ion (X) 


logJV(X) 


[X/H] 


log MX) 


[X/H] 


log MX) 


[X/H] 


Ni 


14.62 ± 0.01 


-1.73 ± 0.07 


> 14.55° 


> -2.08° 


15.25 


± 0.08 


-2.45 


±0.13 


Mgi 


12.40 ±0.15 




12.73 ±0.11 




12.58 


± 0.04 






Mgn 


< 14.69 


< -1.29 


> 14.94* 


> -1.32* 


16.00 


±0.04 


-1.33 


±0.11 


Sin 


15.15 ±0.03 


-0.81 ±0.08 


15.34 ± 0.01 


-0.90 ± 0.05 


15.67 


± 0.03 


-1.64 


±0.10 


Pn 


13.05 ± 0.03 


-0.91 ± 0.08 






13.17 


±0.23 


-2.14 


±0.25 


Sn 


14.66 ± 0.02 


-0.94 ± 0.07 


15.09 ± 0.02 


-0.79 ± 0.05 


15.23 


±0.02 


-1.72 


±0.10 


An 


13.19 ±0.01 


-1.61 ±0.07 






14.42 


±0.02 


-1.73 


±0.10 


Tin 


< 11.85 


< -1.49 


< 12.00 


< -1.62 


12.61 


±0.06 


-2.08 


±0.12 


Cm 


12.87 ± 0.03 


-1.22 ±0.08 


13.13 ±0.06 


-1.24 ±0.08 


13.37 


± 0.01 


-2.07 


± 0.10 


Mnn 


12.35 ± 0.03 


-1.58 ±0.08 






12.84 


±0.03 


-2.44 


±0.10 


Fen 


14.52 ± 0.02 


-1.38 ±0.07 


14.91 ± 0.01 


-1.27 ±0.05 


14.97 


± 0.02 


-2.28 


±0.10 


Nin 


13.43 ± 0.03 


-1.22 ±0.08 


13.82 ± 0.03 


-1.11 ±0.06 


13.70 


± 0.02 


-2.30 


±0.10 


Znn 


12.20 ± 0.07 


-0.87 ±0.10 


12.50 ± 0.03 


-0.85 ± 0.06 


12.79 


±0.02 


-1.63 


±0.10 



" Because of many blends, we take JV(N i) in the central component as a lower limit on the total neutral nitrogen column density. 

* Considering only N(Mg n) in the main component since the other components are either blended with Ly-o- forest absorptions or not detected. 



dLedoux et al.l2003"tlRodrfguez et al.l2006l) . This could be a con- 
seque nce of the relation b etween molecular fraction and metal- 
licity dPetitiean et al.ll2006l) . 

The system toward HE 0027- 1836 is particularly interest- 
ing as it shows absorption from rotational levels J = to 5 (and 
possibly J = 6) in a single well-defined component. The UV ra- 
diation field intensity in which this system is immersed is found 
to be about 20 times the Galactic ambient flux. This could be 
an upper limit as it is estimated assuming that H2 excitation is 
mainly due to UV pumping. However the same value is found 
when estimating the cooling rate in the gas from the C 11* col- 
umn density. 

Thanks to the very high data quality, we observe for the first 
time at high z an increase of the Doppler parameter b from low 
to high H2 rotational levels. To the first approximation, there is 
a linear relation between the kinetic energy of the molecule (as 
given by the Doppler parameter) and the energy of the rotational 
levels. The explanation that this is a direct consequence of the 
formation of H2 onto dust-grains (see lLacour et al.ll2005l) is dif- 
ficult to accommodate given the low formation rate in the DLA 
gas due to low metallicity and small dust-to-gas ratio. 

Detailed PDR models have been constructed to reproduce 
the observations in this system that is special because of the 
low /V(H2,J=2) column density compared to other J level col- 
umn densities. Two components are needed, one with low par- 
ticle density and weak radiation field to reproduce J < 2 levels 
column densities and one with high density and strong radia- 
tion field for J > 2 levels. This combination may appear ad- 
hoc but is unavoidable in the context of PDR models. On the 
other hand, the Doppler broadening and the excitation of H2 
can also be explained by increased temper ature in part of the 
gas due to turbulent dissipation or C-shocks (|Joulain et all 1998; 
ICecchi-Pestellini et al.l 12005; Lac ouret al.l 120051) . as supported 
by the small depletion, the low molecular fraction and the small, 
if any, velocity shift between H2 and C 1 absorptions. 
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Fig. 18. Depletion pattern for the three systems presented in 
this paper. The solid, dashed and dotted lines, represent the typ- 
ical relative abundances observed in, respectively, cold (C) and 
warm (W) gas in the Galac t ic disk and diffuse gas in the Galactic 
halo (H), from lWeltv et all (1 19991) . The dashed-dotted line corre- 
sponds to solar abundances. Filled squares stand for H2-bearing 
components, open squares for other components in the system. 
We show only the components in which Znn is detected. For 
each system, components are shifted in abscissa according to 
their redshift (i.e., the bluest components are on the left). It is ap- 
parent that the Galactic halo depletion pattern represents best the 
observed abundances. There is no significant difference between 
H2-components and the other components although it seems that 
depletion (from Crn, Mnn, Fen, and Nin) is slightly higher in 
the H 2 component of the DLA toward HE 0027- 1836. 
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Fig. 19. Absorption lines from H2 in rotational level J = at 
z abs = 2.40183 toward HE 0027-1836. 
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Fig. 20. Absorption lines from H2 in rotational level J = 1 at 
z abs = 2.40183 toward HE 0027-1836. One can see that damp- 
ing wings are present for H2 L8-0R1. 
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Fig. 21. Absorption lines from H2 in rotational level J = 2 at 
z abs = 2.40183 toward HE 0027-1836. Note the high signal-to- 
noise ratio in the H2 L2-0 R2 and H2 L2-0 P2 panels. 
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Fig. 22. Absorption lines from H2 in rotational level J = 3 at 
z abs = 2.40183 toward HE 0027-1836. 
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Fig. 23. Absorption lines from H2 in rotational level J = 4 at 
z abs = 2.40183 toward HE 0027-1836. The H 2 W0-0Q4 transi- 
tion is actually blended with H2 L7-0 R2 which provides most of 
the absorption. 
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Fig. 24. Absorption lines from H 2 in rotational level J = 5 at 
z abs = 2.40183 toward HE 0027-1836. 
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Fig. 25. Possible detection of H 2 in rotational level J = 6 at 
z a bs = 2.40183 toward HE 0027-1836. A consistent absorption 
feature is seen in the H2L3-OP6 panel, and explains the slight 
asymmetry in the H 2 L4-0 R6 panel, which actually is a blend of 
H 2 L3-0 P2 together with H 2 L4-0 R6. 



